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The possibility has been suggested that high-frequency quasi-periodic oscillations 

O ■ 

■ observed in low-mass X-ray binaries are resonantly excited disk oscillations in de- 
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formed (warped or eccentric) relativistic disks (Kato 2004). In this paper we examine 
this wave excitation process from a viewpoint somewhat different from that of pre- 
vious studies. We study how amplitudes of a set of normal mode oscillations change 



secularly with time by their mutual couplings through disk deformation. As a first 
. step, we consider the case where the number of oscillation modes contributing to the 

resonance coupling is two. The results show that two prograde oscillations interacting 
through disk deformation can grow if their wave energies have opposite signs. 
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1. Introduction 

The origin of high frequency quasi-periodic oscillations (HF QPOs) observed in low mass 
X-ray binaries (LMXBs) is one of challenging subjects to be examined, since its examination 
will clarify the structure of the innermost part of relativistic accretion disks and the spin of 
central sources. One promising possibility is that the QPOs are disk oscillations excited in the 
innermost region of relativistic disks. Particularly, the idea that HF QPOs are disk oscilla- 
tions resonantly excited in deformed (warped or eccentric) disks has been suggested by Kato 
(2004, 2008a,b) by analysical considerations, and studied by Ferreria and Ogilvie (2008) and 
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Oktariani et al. (2010) by numerical calculations. In this model, a disk oscillation (hereafter, 
original oscillation) interacts non-linearly with the disk deformation to produce an oscillation 
(hereafter, intermediate oscillation). The intermediate oscillation is a forced oscillation due to 
the coupling between the original oscillation and the deformation. The intermediate oscillation 
then resonantly responds to the forcing term at a certain radius (Lindblad resonance). After 
having the resonance, the intermediate oscillation interacts non-linearly again with the disk 
deformation to feed back to the original oscillation. Through this feedback process the original 
and intermediate oscillations are excited, if cetain conditions are satisfied. 

Some important consequences have been obtained so far in relation to the origin of 
the instability: i) The wave energies of the original and intermediate oscillations must have 
opposite signs (Kato 2004, 2008a). That is, the instability is a result of wave energy exchange 
between two oscillations with opposite signs of energy through a disk deformation, ii) Since 
the intermediate oscillation resonantly responds to forcing terms resulting from the coupling 
between the disk deformation and the original oscillation, its amplitude is not necessarily small 
compared with that of the original oscillation (see figure 3 of Oktariani et al. 2010). This means 
that the terminology of "original" and "intermediate" oscillations has no particular meaning. 
The two oscillations are equal partners, i.e., the role of the original oscillation mentioned in 
the previous paragraph can be also performed by the intermediate oscillation. Hence, the 
interaction between two oscillations can be schematically sketched as figure 1. This has been 
correctly acknowleged by Ferreira and Ogilvie (2008) (see figure 3 of their paper). 

Based on the above considerations, we develop here a perspective analytical method to 
study the instability. In the previous analytical studies, we considered only the cases where the 
disk deformation is time-independent (i.e., its frequency of the disk deformation, wd, is zero) 
and frequencies of the two oscillations, oj\ and U2, coupling through disk deformation are the 
same, i.e., U\ = uj<i- In the present analyses, the resonant condition is extended to U\ = U2 ±i^d 
with non-zero wd and effects of weak deviation from the resonant condition, u\ = U2 ±wd, on 
growth rate of oscillations are also examined. 

It is noted that the number of oscillations contributing to this resonant process is not 
always limited to two. In order to understand the essence of this instability, however, we 
assume in this paper that only two modes of oscillations contribute to this resonsnt process. 
More general cases will be a subject in the future. 

2. Basic Hydrodynamical Equations and Some Other Relations 

We summarize here basic equations and relations [see Kato (2008a, b) for details] to be 
used to discuss wave couplings through a disk deformation. 
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Fig. 1. Schematic diagram showing two resonant interactions between two disk oscillations (mode 1 and 
mode 2) through couplings with disk deformation. Two oscillations must have opposite signs in their wave 
energy. In this figure, mode 1 is taken to have a negative energy 

2.1. Nonlinear Hydrodynamical Equations 

When we try to quantitatively apply the present excitation process to high-frequency 
QPOs observed in GBHCs (galactic black hole candidates) and LMXBs, effects of general 
relativity should be taken into account. The general relativity is, however, not essential to 
understanding the essence of the instability mechanism. Hence, in this paper, for simplicity, 
formulation is done in the framework of a pseudo-Newtonian potential, using the gravitational 
potential introduced by Paczyhski and Wiita (1980). We adopt a Lagrangian formulation by 
Lynden-Bell and Ostriker (1967). 

The unperturbed disk is in a steady equilibrium state. Over the equilibrium state, 
weakly nonlinear perturbations are superposed. By using a displacement vector, the weakly 
nonlinear hydrodynamical equation describing adiabatic, nonself-gravitating perturbations is 
written as, after lengthly manipulation (Lynden-Bell and Ostriker 1967), 

po^f + 2 P oK ■ V)§ + L(£) = poC(Z,Z), (1) 
where L(£) is a linear Hermitian operator with respect to £ and is 

L{£) = PoK ■ V)(«o • V)£ + po(£ • V)(Wo) + V [(1 - rOpodiv* 

-poV(div^) - V[(£ • V)p ] + (£ ■ V)(Vp ), (2) 

and Po( r ) an d Po( r ) are respectively the density and pressure in the unperturbed state, and Vi 
is the barotropic index specifying the linear part of the relation between Lagrangian variations 
Sp and Sp, i.e., (Sp/p )i inear = r 1 (5p/p )imear- Since the self-gravity of the disk gas has been 
neglected, the gravitational potential, ipo(r), is a given function and has no Eulerian perturba- 
tion. In the above hydrodynamical equations (1) and (2), there is no restriction on the form 
of unperturbed flow u . However, in the followings we assume that the unperturbed flow is a 
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cylindrical rotation alone, i.e., uq = (0,rfi,0), in the cylindrical cordinates (r, tp, z), where the 
origin is at the disk center and the z-axis is in the direction perpendicular to the unperturbed 
disk plane with Q(r) being the angular velocity of disk rotation. 

The right-hand side of wave equation (1) represents the weakly nonlinear terms. No 
detailed expression for C is given here (for detailed expressions, see Kato 2004, 2008a), but 
an important characteristics of C is that we have commutative relations (Kato 2008a) for an 
arbitrary set of ■q 1 , rj 2 , and rj 3 , e.g., 

J po»7i C(r) 2 , r7 3 ) dV = J p ri 1 C{r)z,'n 2 )dV = J Pori 3 C{ri u ri 2 )dV. (3) 

As shown later, the presence of these commutative relations leads to a simple expression of 
instability criterion. We suppose that the presence of these commutative relations is a general 
property of conservative systems beyond the assumption of weak nonlinearity. 

2.2. Orthogonality of Normal Modes 

In preparation for subsequent studies, some orthogonality relations are summarized here. 
Eigen-functions describing linear oscillations in non-deformed disks are denoted by £ a (r,t). 
Here, the subscript a is used to distinct all eigen-functions. Time-dependent part of £ a (r,t) is 
expressed as exp(iu a t), where u a is real. Then, £ a (r,t) satisfies 

-u 2 a p £ a + 2iu aPo (u • V)£ Q + L(U = 0. (4) 

Now, this equation is multiplied by £g(r,t) and integrated over the whole volume, where 
the superscript * denotes the complex conjugate and (3 ^ a. The volume integration of 
Po£,*p( r ,t)£, a (, r it) over the whole volume is written hereafter as (po£p£ a )- Then, we have 

-^<M;0 + 2iu a (pot;(u ■ V)0 + • HU) = 0. (5) 

Similarly, integrating the linear wave equation of £^ over the whole volume after the equation 
being multiplied by £ Q , we have 

-<4(pota&) ~ UMPoUuo ■ V)Q> + (€« • HCp)) = 0. (6) 
Since the operator L is a Hermitian (Lynden-Bel and Ostriker 1967), we have the relation of 

(t a ■ HV P )) = ([L(C)r ■ Q = (HU ■ Q- (7) 

Hence, the difference of the above two equations [eqs. (5) and (6)] gives, when up ^ u a , 

(cu a + UfMpoZJ}) = 2i(po€;(uo • V)$ a ) = -2i(p Q £ a {u ■ V)^). (8) 

In deriving the last equality, we have used an integration by part, assuming that po vanishes 
on the disk surface. 

Different from the case of non-rotating stars, the eigenfunctions of normal modes of disk 
oscillations are not orthogonal in the sense of (po£, a £*p) = 0. In spite of this, however, eigenfunc- 
tions of disk oscillations are orthogonal in many situations. They are classified by azimuthal 
wavenumber, m, node number in the vertical direction, n, and that in the radial direction, 
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£, in addition to the distinction of p- and g-modes [see Kato (2001) or Kato et al. (2008) 
for classification of disk oscillations]. Eigenfunctions with different azimuthal wavenumber are 
obviously orthogonal, i.e., (po$ a £*p) = when m a ^ mp. Even if the azimuthal wavenumbers 
are the same, (po£, a £*p) — when n a ^ rip, if the disk is geometrically thin and isothermal in 
the vertical direction. This comes from the fact that in such disks the ^-dependence of eigen- 
functions with n node(s) in the z-direction (n is zero or a positive integer) is described by the 
Hermite polynomials "H n as 

£r, ^OcUniz/H), ^txHn-^z/H), (9) 

(Okazaki et al. 1987), where the subscripts r, ip, and z represent the cylindrical coordinates (r, 
if, z) whose origin is at the disk center and the z-axis is perpendicular to the disk plane. Here, 
T-L n is the Hermite polynomial of argument zj H, H being the half-thickness of the disk. Thus, 
the eigen-functions classified by m and n are orthogonal. 1 In summary, we have 

(10) 

where 5 a) b is the Kronecker delta, i.e., it is unity when a — b, while zero when a ^ b. 

Orthogonality of (po£, a £*p) does not hold in the case where m a = mp and n a = rip. Even 
in these cases, however, (po$, a £*p) will be close to zero for £ a ^ £p, if our interest is on short 
wavelength oscillations in the radial direction, since the radial dependence of eigenfunctions is 
close to sinusoidal in such cases. 

3. Couplings of Two Oscillations through Disk Deformation 

Let us consider the case where two oscillation modes, £i(r,t) and £ 2 ( r >t)i resonantly 
couple through disk deformation ^D( r ^)- Through the coupling term C(£,£) [see eq.(l)] many 
other modes than £ x and £ 2 appear, and their amplitudes as well as those of ^ l and £ 2 become 
time dependent. Now, we assume that normal modes of oscillations form a complete set, and 
expand the resulting oscillations, £(r,t), including the disk deformation, £ D (r,t), in the form 

^r,t) = A 1 (t)^ l (r,t)+A 2 m 2 (r,t) + A D ^ D (r,t)+Y,A a (t)Ur,t). (11) 

a 

Since our main concern is on the modes 1 and 2, £ x and £ 2 are distinguished from other eigen- 
functions and the subscript a is hereafter used only to denote other eigenfunctions than £ x 
and £ 2 - O ur purpose here is to derive equations describing a secular time evolution of A\ and 
A 2 . The disk deformation £v(r,t) is assumed to have a much larger amplitude than other 
oscillations and its time variation during the coupling processes is neglected, i.e., A-q = const. 

We now express the eigen-frequencies associated with £ 1; £ 2 , £ D , and £ a by cux, u 2 , uj^, 
and u a , respectively, i.e., £x(r,i) =exp(ia>it)£ 1 (r) etc. Then, substitution of equation (11) into 
equation (1) leads to 

1 In vertically polytropic disks, orthogonality holds by using the Gegenbauer polynomials (Silbergleit et al. 
2001). 
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i=l,2 z 



(u - V)]^ 1 + 2p ^[^2 



/7 4 

(«o • V)]£ a + ]T2p ^[zu; Q + (u • V)]£ c 



(12) 



where terms of d 2 A 1 /dt 2 , d 2 A 2 /dt 2 , and d 2 A a /dt 2 have been neglected, since we are interested 
in slow secular evolutions of A's. On the right-hand side of equation (12), the coupling terms 
that are not related to the disk deformation are neglected. 2 

Now, we define the wave energy E\ of normal mode of oscillation £ x by 
1 



E 1 



-uji 



(13) 



(Kato 2001, 2008a). To use later, the wave energy E 2 of the normal mode £ 2 is a l so defined by 



E 2 = -0J2 



Furthermore, we introduce the following quantities: 
1 



W 12 = 
W 12 * 



W lM = -({p ? 1 -C{£ a ,&] 



(14) 

(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



To proceed further, the time and azimuthal dependences of normal mode oscillations are 
written explicitly as 

£fc( r ^) = exp[i(u k t - m k p)]£ k (k = 1, 2, a, D). (21) 

To avoid writing down similar relations repeatedly, we have introduced the subscript k, which 
represents all oscillation modes, i.e., k denotes 1, 2, a and D. Here, we take all m k to be zero 
or positive intergers, while u k is not always positive. If u k < 0, the oscillation is retrograde. 



In writing down the right-hand side of equation (12), we have used the following relation: 

K(i4)K(J3) = ^[AB + AB*], 
where A and B arc complex variables. 
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It is noted here that by using equation (21), we can express the wave energy of the 
normal mode oscillation in an instructive form. Since the r- and </?- components of £ 1; say £i r 
and £1^, are related in a geometrically thin disks by (e.g., Kato 2004) 

i(wi-mfii)^ + 2fi^i r ~0, (22) 

we have 

E 1 ~ ^({ux -rmtypoitUir + (23) 



2 

This shows that the sign of wave energy is determined by the sign of uj\ — mitt in the region 
where the wave exists predominantly (e.g., Kato 2001). For example, a prograde (aj\ > 0) wave 
inside the corotation resonance has a negativce energy, while a prograde wave outside it has a 
positive energy. 

In previous papers we mainly considered the case of ui\ = UJ2 with wd = 0. In this paper, 
we extend our analyses to more general cases of resonance: 

uj\ ~ 0J2 ± wd, (24) 

where cud is not necessary to be small. We introduce A + and A_ defined by 

A + = uji — oj 2 — c^d an d A_ = U\ — u 2 + wd- (25) 

In the resonance of u\ ~ U2 + wd, A + is small, but A_ is not small unless wd is small. In the 
resonance of u)\ ~ u 2 — ud, on the other hand, A_ is small, but A + is not always so unless wd 
is small. It is noted that the resonant condition concerning azimuthal wavenumber is 

mi = m 2 ±mu, (26) 

where mi and m 2 are zero or positive integers, while mo is a positive integer, since we focus 
our attention only on the case where the disk deformation is non-axisymmetric, i.e., mo 7^ 0. 

After these preparations, we integrate equation (12) over the whole volume after multi- 
plying £*(r,t). Then, the term with dAi/dt becomes i(4Ei/ui)(dAi/dt), while the term with 
dA 2 /dt vanishes since m 2 7^ mi by definition. Concerning the term with dA a /dt, some more 
consideration is necessary. If uj a ^u>i, by using equation (8) we can reduce the integration to 
d A 

*£^K-^i)<MtO- (27) 

In the case where m a = mi im^, (Po£i£a) vanishes since 7^ 0. On the other hand, when 
m a 7^ mi ±771d, A a does not appear in the coupling term. That is, A\ and A a have no nonlinear 
couplings and thus we can take as A a = when we consider time evolution of A 1 . In the case 
of ui a = oji, the last term on the left-hand side of equation (12) does not lead to equation (27). 
Even in this case, by the same arguments as the above we can neglect the term. Considering 
these situations we have 

dAi 4E 1 
i— 

dt u>i 



= A 1 (A B W n + A* D W n *) + A 2 (A D W 12 + A* D W 12 *) + ]T A a (A D W la + A* D W la *). (28) 

a 

Among various coupling terms on the right-hand side of equation (28), the first two 
terms with Wu and Wu* can be neglected if we consider non-axisymmetric disk deformations 
such as warp or eccentric deformation, since the terms inside ( ) in equations (15) and (16) are 
proportional to exp(— imutp) and exp(imB<p), respectively, and their angular averages vanish. 
The last two terms with W\ a and W\ a * on the right-hand side of equation (28) are also neglected 
hereafter by the following reasons. The terms inside of ( ) of equations (19) and (20) are 
proportional to exp[i(— uj\ + u a + wo)t] and exp[z(— u\ + uj a — cud)£], respectively. In general, 
they rapidly vary with time, since no resonant condition is assumed among uj\, u a , and wd- 
Hence, if short timescale variations are averaged over, 3 the averaged quantities are small and 
can be neglected. 4 In summary, the coupling terms remained are the middle two terms of 
equation (28) with Wi 2 and Wi 2 *. 

Based on the above preparations, we reduce equation (28) to 



^—-jr = A 2 A D W 12 exp(-iA + t)5 mum2+mu 

+A 2 A* D W 12 *exp(-iA_t)5 (29) 



where the symbol 5 a ,b means that it is unity when a — b, but zero when a^b. Here, from W\ 2 
and W\2* the time and azimuthally dependent parts are separated as 

W 12 = W 12 exp(-iA + t)5 mum2+mD (30) 

Wi 2 * = H / i2*exp(-iA_t)5 miim2 _ mD . (31) 

A physical meaning of equation (29) is as follows. The imaginary part of (ui/2)Wi 2 , for example, 
is the rate of work done on mode 1 (when mode 2 and the deformation have unit amplitudes) 
through the coupling of mi = mi +mn (Kato 2008a). Hence, in a rough sense, equation (29) 
represents the fact that the growth rate of mode 1 is given by energy flux F-±[= A 2 Ad(ui/2)Wi 2 } 
to mode 1 as 

^i = ^. (32) 
dt 2E 1 K ' 

The next subject is to derive an equation describing the time evolution of A 2 by a similar 

procedure as the above. That is, equation (12) is multiplied by £2 an d integrated over the whole 

volume. Then, as the equation corresponding to equation(28), we have 



We are interested in solutions where all A's vary slowly with time 

In some cases, however, some of u> a 's is/are close to u)\ ±wrj. For example, all trapped axi-symmctric g-modc 
oscillations have frequencies close to At ma x (the maximum of the epicyclic frequency), when their azimuthal 
wavenumber is zero. Then, W\ a or W\ a „ have no rapid time variation and cannot be neglected by the above 
argument of time average when m a satisfies the relation of m a = m\ ±tod- Then, the terms with W\ a 
or W\ alf also contribute to resonant couplings. Such cases are outside of our present concern. See related 
discussions in the final section. 
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dA 2 4E 2 

i— 

at u 2 

= A 1 (A D W 21 + A^W 2U ) + A 2 (A D W 22 + A* D W 22 *) + Y. A MvW 2a + A* D W 2a *). (33) 

a 

Here, not all of the expressions for W 2 i, W 2 u, W 22 , W 22 *, W 2a , and W 2a * are given, since they 
are the same as those of W\\, W\u, W± 2 , Wi 2 *, Wi a , and Wi a *, respectively, except that £* in 
the latters are replaced now by £ 2 . As examples, we give only W 2 \ and W 2 i* as 

W21 = \ ((p £ ■ C(€ 1) fo))+(pbG ■ C(€d,*i))), (34) 

^21* = \ ((po^ ■ Cfa^+fpoS ■ C(&,tS)) ■ (35) 

By the same argument used in reducing equation (28) to equation (29), we simplify 
equation (33). That is, in the present case the coupling terms that remain are those with W21 
and W 2 u, and we have finally 

4i— = AiA D W2iexp(z'A_t)5 m2imi+mD 
UJ 2 Cut 

+A 1 A* D W 2U exp(iA + t)5 m2 , mi _ mo , (36) 

where W 2 i and W\ 2 * are the time and azimuthal dependent parts of W 2 \ and W 2 \*: 

W 21 = W 21 exp(iA^t)S m2 (37) 

W 2U = W 2 i*exp(iA + t)6 m2tmi - mD . (38) 

The main results obtained in this section are equations (29) and (36). 

Here, it is of importance to note that we have the following identical relations: 

W 2l = (Wia.)*, (39) 

W21* = (W 12 )*, (40) 

where the superscript * means the complex conjugate. These relation come from the fact that 
for arbitrary functions, rj 1 , r] 2 , and rj 3 , their order in (por) l ■ C(r) 2 ,r) 3 )) can be arbitrary changed 
[see equation (3)] (Kato 2008a). 

4. Growth Rate of Resonant Oscillations 

By solving the set of equations (29) and (36), we examine how amplitudes of A\ and A 2 
evolve with time. We consider two cases of m 2 = vii\ +mj) and m 2 = mi — Wd, separately. 

4-1. Case of m 2 = mi + rriD 

In this case the set of equations of A± and A 2 are, from equations (29) and (36), 

A .ExdA± = A2A ^w 12 *eM-iA_t), (41) 
u)\ Cut 
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4 .E2 dM = AlAD w 21 exp(iA_t). (42) 

id 2 dt 

By introducing a new variable A\ defined by 

A 1 = Aiexp(zA_t), (43) 

we can reduce the above set of equations to 

4 .£i dM + 4 ^i A _ ii = 4,4^ (44) 
Wi at tui 

4*^^ = ^ D W 21 . (45) 
u 2 at 

Hence, by taking A\ and A 2 to be proportional to exp(iat), we obtain an equation describing 
cr as 

where equation(39) is used. 

In the limit of an exact resonance of A_ = 0, the instability condition (cr 2 < 0) is found to 
be (idi/Ei)(u 2 /E 2 ) < 0. The meaning of this condition is discussed later. If the frequencies of 
two oscillations deviate from the resonant condition of u\ = U2 — wd, the growth rate decreases. 
This can be shown from equation (46). That is, the condition of growth is 

A 2 _ + ^|A D | 2 |W 12 *| 2 <0, (47) 

and the growth rate tends to zero, as A 2 , increases from zero. If A 2 increases beyond a certain 
limit the left-hand side of inequality (47) becomes positive, and a is no longer complex. That 
is, the amplitude of oscillations are modulated with time, but there is no secular increase of 
them. 

4-2. Case of m 2 = m\ — m^. 

In the present case, from equations (29) and (36), we have 

4i — = A 2 A D W 12 exp(-iA+t), (48) 
oj 1 at 

u E*dM = w 2U exp{iA + t). (49) 

id2 at 

A new variable Ai is introduced here by 

Ax = A^xpiiA+t). (50) 
Then, the set of equations (48) and (49) are reduced to a set of equations of A\ and A 2 as 

dM + m ^ = (51) 

id\ dt id\ 

^E 2 dM =A ^ 

id 2 dt 
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Hence, by taking A\ and A 2 to be proportional to exp(zat), we have 

^"^"S^ 1 ^ 1 ^ 0, (53) 

where we have used equation (40). 

Two oscillations certianly grow again at the limit of exact resonance of A + = (i.e., 
U\ = u 2 +wd), if (^i/Ei)(u 2 / E 2 ) < 0. Even if the resonance is not exact, they grow if is 
small enough so that 

is satisfied. 

4-3. A Relation between A\ and A 2 

Finally, it is useful to derive an instructive relation between A± and A 2 . Let us first 
consider the case of m 2 = m\ +wd. Let us multiply A\ to equation (41) and also A\ to the 
complex conjugate of equation (41). Then, summing these two equations we have an equation 
describing time evolution of |>li| 2 . Similarly, from equation (42), we can derive an equation 
describing time evolution of \A 2 \^ . Summing these two equations, we have finally 



d_ 
dt 



— A 1 2 + — A 2 2 

U)i u 2 



0, (55) 



where we have used W 2l = W\ 2 *. The same equation can be derived from euqations (48) and 
(49) in the case of m 2 =m\ — m D . To derive the equaition, W 2 u = W{ 2 has been used. 

Equation (55) obviously shows that (yj\f Ei)(u 2 / E 2 ) < is necessary for growth of os- 
cillations. In the case of (ui/Ei)(u 2 /E 2 ) > 0, on the other hand, amplitudes of A 1 and A 2 
are limited, although the relative amplitude of both oscillations may change with time by 
interaction through disk deformation. 

4-4- Summary of Resonant Instability Condition 

The results in the previous subsections show that when resonant conditions of u>\ = 
u 2 and mi = m 2 ± are satisfied among two oscillations characterized by (ui, mi) and 
(u 2 , m 2 ) and disk deformation characterized by (cud, tid), the two oscillations are resonantly 
excited if (ui/Ei)(u 2 /E 2 ) < is realized. A deviation from the condition of u>i = u> 2 ± u>d 
decreases the growth rate, but oscillations grow as long as the deviation is smaller than a 
critical value. 

In the case where both of U\ and u 2 are positive (i.e., both oscillations are prograde), 
the above instability condition is EiE 2 < 0. This is the result suggested by Kato (2004, 2008a, 
b) by a different approach. 

In the case where wd is larger than u 2 (> 0), a resonant condition, cui = u 2 -ud, is 
satisfied for u>i < (i.e., retrograde wave). In this case the wave energy E\ is positive [see 
equation (23)] and the instability condition is reduced to E 2 /u 2 > 0, i.e., E 2 > 0. That is, when 
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uJiU2 < 0, the condition of resonant instability is E\Ei > 0. It is noted that in the case where 
both of u\ and LO2 are negative, both E\ and E2 are positive, so that the condition of resonant 
instability, (ui/ Ei)(u2/ E 2 ) < 0, cannot be satisfied. 

Among three case of i) U\ > and 0J2 > 0, ii) Ui < and u>2 > 0, and iii) U\ < and 
< 0, the interesting case in the practical sense is the first one, which is discussed in the next 
section. 

5. Discussion 

First, let us describe, in terms of the present formulation, the g- and p-modes resonant 
instability that was numerically studied by Ferreira and Ogilvie (2008) and Oktariani et al. 
(2010). They considered the resonant interaction, through a standing warp (w D = 0), between 
i) the axisymmetric g-mode oscillation whose £ r has one node in the vertical direction and ii) 
the one-armed p-mode oscillation whose £ r has no node in the vertical direction. That is, the 
set of (u, m, n) is (~ K max , 0, 1) for the g-mode oscillation, and (~ ft max , 1, 0) for the p-mode 
one, where /t max is the maximum of the (radial) epicyclic frequency. The warp is taken to be 
(0, 1, 1). In this case, the resonant conditions, u>i ~ u 2 > (i.e., w D = 0) and m 2 = mi +mD, 
are satisfied. Hence, if EiE 2 < 0, these modes are excited simultaneously. This condition of 
E1E2 < is really satisfied, since the axisymmetric g-mode oscillation has a positive energy, 
while the p-mode oscillation trapped between the inner edge of the disk and the barrier resulting 
from the boundary between the propagation and evanescent regions has a negative energy. The 
cause of excitation is a positive energy flow from a negative energy oscillation (p-mode) to 
a positive energy oscillation (g-mode). By this energy flow both oscillations grow. The disk 
deformation is a catalyzer of this energy flow. 

In the above argument the disk deformation is assumed to be a warp. Instead of a warp, 
we can consider a c-mode oscillation as one of other possible disk deformations. In this case 
the set of (u, m, n) of the disk deformation is (wd, 1, 1), where ujd is the frequency of c-mode 
oscillation and wd <S K max , unless the spin of the central source is high (Silbergleit et al. 2001). 
The resonant condition in this case is u\ = U2 ± ud, not u\ =002- 

There are some limitations of a direct comparison of the present analytical results to the 
numerical ones by Ferreira and Ogilvie (2008). In the present analyses only two normal modes 
of oscillations are assumed to contribute to the resonance to understand the essence of the 
resonant instability. In the realistic case considered numerically by Ferreira and Ogilvie (2008), 
however, more than two normal modes of oscillations may contribute to the resonance. In their 
case one of resonant oscillations is an axisymmetric g-mode. As mentioned in footnote 4, eigen- 
frequencies of the axisymmetric (m = 0) g-mode oscillations with different I (with n — 1) are all 
close to K max . Hence, the g-mode oscillations that satisfy the resonant conditions, u± =U2 + ^d 
and rrii — m 2 + ^d ; m ay not be only one, and overtones of g-mode oscillations with nodes in 
the radial direction may also contribute partially to the resonance. If this is the case, some 
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coupling terms with other A than A\ and A2 appear in equations describing time evolution of 
A\ and A2. In the present analytical formulation such situations are not considered. To extend 
our analyses to such cases, we must derive equations describing the time evolution of other 
A's than A\ and A2, and these equations should be solved simultaneously with the equations 
describing time evolution of A\ and A2. We think that the essence of the instability mechanism 
is already presented in the case where couplings occur only between two oscillations. However, 
since such extention of our formulation is formally simple and there may be some subsidiary 
modifications of instability criterion, such extension should be done in the near future. 

In the present formulation, some important theoretical problems remain to be clarified. 
One of them is whether the set of normal modes of oscillations form a complete set. If not, it is 
uncertain whether the oscillations realized on the disk can be expressed in the form of equation 
(11). 

In the case where wd is larger than U2, the resonant frequency uj\ which satisfies the 
condition, oj\ = cj 2 — wd, is negative. In this case of uj\ < 0, E\ is positive [see equation (23)]. 
Hence, the instability condition in this case is E\E^ > 0, as mentioned in the last section. All 
the g-mode and p-mode oscillations with negative frequency are, however, not trapped in the 
inner region of disks. They propagate away far outside unless the disks are truncated, 5 That 
is, their frequencies are continuous and will not be directly related to the QPO phenomena in 
disks. 

It is important to note here that the innermost region of relativistic disks is a place 
where the present excitation mechanism works most efficiently. For the mechanism to work, 
two normal modes of oscillations with opposite signs of wave energy must coexist in a common 
region of disks. In general, for m^O, positive energy oscillations propagate in the region outside 
the radius of corotation resonance, while negative energy oscillations inside the resonance. 
Hence, there is a tendency that the propagation region of oscillations with opposite signs of 
wave energy are spatially separated, unless different types of oscillation modes are considered. 
If the propagation region of oscillations with opposite signs of wave energy are separated, the 
coupling efficiency between the two oscillations are weak, i.e., W's in equations (29) and (36) 
are small, and practically there is no growth of oscillations. For oscillations with opposite signs 
of wave energy to coexist in a common region, effects of general relativity is important. It is 
noted that in such cases there is a tendency that the Lindblad resonance of an oscillation (e.g., 
p-mode) occurs in the propagation region of the other oscillation (e.g., g-mode). This helps to 
increase the growth rate since the coupling terms become large, as mentioned before. 

Finally, as an application of the present resonant excitation process of oscillations, we 
briefly note Kato's model of high-frequency twin QPOs (e.g., Kato and Fukue 2006). In this 
model the lower-frequency QPO of the twin is related to the set of g-mode and p-mode os- 

5 This can be shown by examining local dispersion relation of these oscillations. The vertical p-mode oscilla- 
tions that are trapped in finite region also have u)\ > 0. 
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cillations considered numerically by Ferreira and Ogilvie (2008) and Oktariani et al. (2010). 
The higher-frequency QPO of the twin is considered to be the set of one-armed g-mode and 
two-armed p-mode oscillations [see for details figures 1 to 3 of Kato (2008a)]. Here, the g-mode 
oscillation has a positive energy and the p-mode one has a negative energy, and thus this set 
of oscillations can also grow by the resonant coupling. 

In his model, correlated time variation of twin QPOs observed in neutron-star low-mass 
X-ray binaries are described by assuming time variation of uu- One of problems of this model, 
however, is that masses of neutron stars required to take into account observations are rather 
high, i.e., for example, 2.4 M & for Sco X-l and 4U 1636-53. See Lin et al. (2010) for detailed 
comparison of various QPO models with observations. In black-hole low-mass X-ray binaries, 
the twin QPOs have no frequency change with frequency ratio of 3 : 2. In Kato's model this 
can be described by assuming wd = 0. The spin parameters a*'s estimated by this model for 
black hole sources with measured masses are around a* ~ 0.4 (e.g., Kato and Fukue 2006). 
However, the spin parameters estimated by comparing model continuum X-ray spectra with 
observations are generally higher than the above, say around a* ~ 0.8 or more (e.g., Narayan 
et at. 2008). 
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